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ABSTRACT 

Weak lensing calculations are often made under the assumption of the Born approximation, 
where the ray path is approximated as a straight radial line. In addition, lens-lens couplings 
where there are several deflections along the light ray are often neglected. We examine the 
effect of dropping the Born approximation and taking lens-lens couplings into account, for 
weak lensing effects up to second order (cosmic flexion), by making a perturbative expansion 
in the light path. We present a diagrammatic representation of the resulting corrections to the 
lensing effects. The flexion signal, which measures the derivative of the density field, acquires 
correction terms proportional to the squared gravitational shear; we also find that by dropping 
the Born approximation, two further degrees of freedom of the lensing distortion can be ex- 
cited (the twist components), in addition to the four standard flexion components. We derive 
angular power spectra of the flexion and twist, with and without the Born-approximation and 
lens-lens couplings and confirm that the Born approximation is an excellent approximation 
for weak cosmic flexions, except at very small scales. 
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1 INTRODUCTION 

Wea k gravitational flexions are third order distortions of the images of distant galaxies due to gravi tational lensing by the large-scale structure 
jGold berg & Natarajan 2002; Goldberg & Bacon 2005; Bacon et al. 2006; Schneider & Ei 2008[). Flexions have the potenti al to probe the 
profile jLaskv & Fluke 2009; Velander et al. 2010; Fluke & Laskv 2011), substructure jRowe et alJl2007l;lBacon et aljlioioh and ellipticity 
of dark matter haloes ( .Er & Schneidern201Q ; .Hawken & Bridle,.2009,) . improve mass reconstructions i Er et alfcOlOh . or complement cosmic 
shear studies as they are li kely to be less prone to intrinsic alignment effects (for a review of those, see S chaf eil 2009|) . Thev can be measured by 
a shapelet decomposition (Massev et al. 2005; Massev & Refregier 2005 ; Massev et al. 2007; Goldberg & Leonard 2007) or by the octupole 
moments (Irwin & Shmakova 2005, 2006; Irwin et al. 2007; Okura et al. 2007, 2008; Umetsu et al. 2008; Okura & Futamase 2009J of the 

brightness distribution of the lensed galaxies. 

A common approximation in the theory of weak lensing (for reviews, see ISchneideretai]|l992l ;l Bartelmann &Schneideil 1200 ll; 
lBartelmanrj2010l) . in particular in weak cosmic shear, is the Born-approximation in which one integrates the gravitational tidal field (weighted 
suitably for lensing) along a fiducial straight ray instead of the actual photon geodesic. The Jacobian matrix, which describes the image map- 
ping to first order, is symmetric in this approximation. It becomes asymmetric, howe ver, if corrections from the perturbed geodes ic are 
taken into account, and/or if multiple deflections along the line of sight are considered jBernardeau et all 1997; Sch neider et al.ll 19981). The 
infl uence of perturbed geod esies on the cosmic shear signal was worked out in detail by ICoorav & Hiil ( |2002) . .Shapiro & Cooravl ilOOdt) 
and lKrause & Hiratal j200^ , who in additio n treat the complication t hat the actual observable in lensing is the reduced shear g = yl{l - k) 
instead of the gravitational shear y. Recentlv. lBemardeau et alj j2010h pursue a different approach using the Sachs-equation, which describes 
the infinitesimal changes in the cross section of a light bundle during propagation; they derive corrections to the power spectra due to the 
perturbed geodesic and identify further relativistic effects in weak lensing image distortions. These studies ag ree that correct ions due to 
perturbed geodesies are minor, but that vortical p atterns in the shear field can be excited, referred to as S-modes ('Stebbins"l996^. Addition- 
ally, rotation of galaxy images can occur (see e.g. lHirata & Seljak 2003; Pen & Mao 2006). The validity of perturbative techniques has been 
checked against ray-tracing on N-body simulations of the cosmic density field, and found to be well approximated even on nonlinear scales 
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jjain et al.l26ool : l\in Waerbeke et al .l200ll:lHilbert et al.l2008l:lTevssier et alj2009l) . Statistics of Bom-approximated shiear and flexion fields 
have been investigated in detail in lMunshietal.l ( l2010t) and lMunshi et al. 

In this paper, we aim to generalise these results by examining the elTect of dropping the Born approximation, and including multiple 
deflections, at the order of gravitational flexions, which is significantly more complex. Flexions describe the components of the tensor 
didjdktf/ i.e. the third angular derivatives of the lensing potential if/; the Born approximation forces this tensor to be symmetric, due to the 
interchangability of any pair of partial derivatives, giving 4 degrees of freedom. The remaining 2 degrees of freedom, coined twist, can be 
present as a systematic or can be physically excited by dropping the Bom-approximation, which as we will see destroys the symmetry of 
the te nsor, didjdkijj djdidi,i/f. The resulting full set of six flexion-related image distortion modes has been described bv lBacon & Schafe3 
( 12009 ). We aim to provide corrections to the spectra of all flexion components due to the perturbed photon geodesic and to quantify the 
correlation properties of the newly emerging flexion degrees of freedom. While carrying out the computation, we will point out the analogies 
to gravitational shear. We work in a flat-sky approximation, because gravitational flexions are a small-scale phenomenon. 

The paper is organized as follows. After a compilation of the relevant formulae concerning cosmology and structure formation in Sect.|2] 
we introduce gravitational flexions and compute the correction terms due to relaxation of the Bom-approximation in Sect. [3] We also present 
a diagrammatic representation of the distortion corrections. Angular power spectra of the flexion quantities with their corrections are shown 
in Sect.|4] and the main results are summarised in Sect.|5] 



2 COSMOLOGY 

In this section we introduce the quantities we require in order to calculate cosmological weak lensing statistics. We require a description of 
both the cosmological background and the matter fluctuations in the Universe, and we need to be able to describe how these fluctuations 
grow with time. 



2.1 Dark energy cosmologies 

First, we need to describe the background expansion of the Universe. In spatially flat dark energy cosmologies with matter density parameter 
n„, the Hubble function H(a) = dlna/dt is given by 

HHa)_^[l^ 1 - Q,„ 

with the dark energy equation of state parameter w. The value w = - 1 corresponds to a cosmological constant A. The relation between 
comoving distance x and scale factor a is given by 

r' da 

in units of the Hubble distance = c/Hq. 

The cosmological model used throughout is a spatially flat ACDM cosmology with Gaussian adiabatic initial perturbations in the cold 
dark matter density field. The specific parameter choices are fl,,, = 0.25, = 1, erg = 0.8, O4 = 0.04, w = -1 and Ho = 100 h km/s/Mpc, 
with h = 0.72. For simplicity, we assume that the lens ed galaxies reside at a redshift of Zs = 0.9, corresponding to the projec ted median of 
the redshift distribution of the EUCLID galaxy sample jRefregier & the DUNE coUaborationll2008l : lAmara & Refregieill2007h . 



2.2 Matter power spectrum 

We also need to describe the density fluctuations in the Universe. The linear cold dark matter (CDM) density power spectrum P(k) describes 
the fluctuation amplitude of a Gaussian statistically homogeneous and isotropic density field 6, {6(k)6(k')) = {InYSoik + k')P(k), and is 
given by the ansatz 

P{k)ock"'T\k), (3) 
with the transfer function T{k). In low-£l,„ cosmologies T(k) can be approximated by ( lBardeenetai]|l986h . 

T{q) = ^"^^2^34^^^^ ^ ^■^'^^ ^ ^^^-^^^^ ^ ^ (6.71#)"' , (4) 

where the wave vector k = qT is rescaled with the shape parameter F dSugivarn^l 1 995b describing corrections due to the baryon density Q.1,, 

F = Q,„AexpJ-ni,|l + ^jj. (5) 

This fit to T(k) is sufficient here for illustrative purposes. The spectrum P(k) is normalised to the standard deviation of matter fluctuations erg 
on a scale = 8 Mpc/A, 

If C 2 k^ 

a-l = — I dkk'^P(k)W\kR) = I dlnk[A(k)W(kR)) with A^{k)=—P{k) (6) 

where W i s a Fourier transformed spher ical top hat filter function, W{x) = 3 ji(x)/x, and j((x) is the spherical Bessel function of the first kind 
of order £ jAbramowitz & Stegurj|l972h . 
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Figure 1. Image distortions induced in gravitational lensing: a constant potential does not affect the image of a galaxy (first panel); a potential varying linearly 
with position displaces an image (second panel); a potential with second derivatives shears an image (third panel); and non-vanishing third derivatives of the 
potential bend an image (fourth panel), described by the flexions. Hence, the flexion is related to a quadratic variation of the displacement, or a linear variation 
of the shear across the galaxy image. 

2.3 Structure growth 

Finally in this section, we need to describe the growth of structure as time progresses. Linear homogeneous grow th of the density field , 
6{x,d) = D+(a)d(x, a - 1), is d escribed by the gr owth function D+(a), which is the solution to the growth equation jTumer & Whitd[l997l : 
IWang & Steinhardtll 19981 : iLinder & Jen]dnil2003l) , 

-^Z)^(fl) + -(3 + ^D^(a) = -^fi„(a)D+(fl), (7) 

aa-^ a\ dlna / da 2a^ 

whose solution is D+(a) = a if Q.„, = 1 and remains close to this solution in dark energy cosmologies. Consequently, the spectrum grows 
according to P{k,a) = D\{a)P{k). For the purpose of this paper we restrict ourselves to this linear growth regime, despite flexion having 
much higher signal-to-noise on smaller scales. This is due to the computational complexity of the calculations below; however, we will see 
that the correction features can be easily understood in this regime and quantitatively extrapolated correctly to the nonlinear regime. 

3 GRAVITATIONAL LENSING 

In this section we will show how the lensing quantities up to third order can be calculated, including lens-lens couplings and dropping the 
Born approximation. We will also introduce a scheme of diagrams that can describe the various lensing correction terms. 



3.1 Weak lensing and the light propagation equations 

Weak gravitational lensing occurs when the shapes of distant galaxies are distorted as their light is deflected by intervening gravitational 
potentials. The light deflection translates the image of a galaxy to a new position on the celestial sphere. If the deflection varies linearly 
across the galaxy image, it appears sheared. If furthermore the deflection shows an appreciable quadratic variation across the image (i.e. 
the shear varies linearly with position), one can observe a bending of the image, which is refered to as gravitational flexion. Fig.[T]gives an 

impression of the image distortions in weak gravitational lensing. 

In order to describe these phenomena in detail, we use the common formalism telandfordetd] Il99ll ; ISeitz et alj 1 19941 ; 
ISeitz & Schneidei|[r994 ; .Kaiser. 1996; Jain & Seljak 19971 ; [Schneider et al.lll998h for describing light propagation in an inhomogeneous 
universe: x{0,x) is the transverse comoving separation between a light ray and a fiducial ray which intersect at the observer at an angle 6. 
This separation increases linearly with comoving distance and acquires corrections due to light deflection in gravitational potentials. 



(8) 



x(6,x)=xd-- dy (A'-y)[v^o(A:(0,y),y)-v^ci)(''>Cr')]=Te-2 dy (x-x')['^MxiO,x'),x')-'^^v'^°Hx')], 

where we have defined the dimensionless potential <p = for brevity. Without lensing, a source at comoving distance separated by x 
from the fiducial ray, will be observed at an angle (3 = x/x- With lensing, this angle becomes 



/i,(0,x) = e,-2 Ax 



,x-x 



X 



U,{x{e,x'\x')-vTix') 



(9) 



where the index on the gravitational potential denotes the derivative perpendicular to the line of sight relative to the component / of the 
physical coordinate. The shear and convergence of galaxy images is related to the Jacobian matrix = d/3/d6 of the mapping, which can 
be obtained by differentiation of the lensed source position (3 with respect to the true celestial coordinate 6 of the source: 

d/3, 



^„(B,X) = 



^ S.j - 2 



dx' Or -x')-v,kix{e,x'lx'mji0,x'l 

X 



(10) 



with summation over repeated indices. Note that this will be quite difficult to solve in general, as it is an implicit relation, with ^ appearing 
on both sides. In order to examine flexion-related image shape changes, we require the angular derivative of the Jacobian ^ which defines 
the rank-3 tensor S: 
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Jo 



X 



iPunAx{e,x\x')^uiO,x)^mk{e,x') + -Vun{x{9,x\x')S,,,jk{e,x') 

X 



(11) 



It is worth noting that is always symmetric under exchange of the last two indices, S,ji- = which follows straightforwardly from its 
definition as the partial derivative of JTI: 



S, 



•ikj- 



Due to this symmetry there are in total 6 independent entries in Sijk'- Sooo, Soio = 2?ooi, Sou, 2?ioo, Sioi = Suo and Sm. 



(12) 



3.2 Perturbative solution to the optical equations 

The equations for the lensing quantities above can be solved perturbatively: starting with the solutions for the unlensed case (i.e. a straight 
line, which corresponds to the Born approximation) one recoveres iterative corrections to the geodesic x(9,x) and hence for the image 
position I3{6,x), for the Jacobian y{(0,x) and finally for the derivative S(0,x)- The perturbative expansion is controlled by considering 
corrections to the geodesic up to 0{(p) in the gravitational potential. Consequently, the Jacobian acquires corrections of order 0(<p^) and the 
derivative is corrected by terms of order 0(<p^). 

The photon geodesic is approximated by x = jc*"' + jc*'' + 0(tp-), with the zeroth order solution x^^'^OjX) being a straight line x'^''''(9,x) = 
xd, and a first order correction x'-'\0,x) which is linear in the gravitational potential, as seen in equation l[8j: 

x^'H0,x) = -2 r dx' (x-x')['^Mx°(e,x'),x') - V^(^<'"Ct')] . (13) 

Jo 

The Jacobian matrix is expanded to second order in the gravitational potential, ^ = ,51*°' + ,51*'' + ,51*^' + 0((p^), and has the unit-matrix as the 
zeroth-order solution J{fj^ = 6ij which corresponds to the unlensed case, i.e. the identity mapping. The first order correction J?I*j' measures 
the line of sight integrated tidal fields (fjj (the indices on the gravitational potential denote a partial derivative): 

Ji'^(e,x) = -2 r dy ix -x')-v.M\e,x'\x'\ (i4) 

-' Jo X 

and the second order correction Jif^ comprises terms caused by lens-lens coupling (oc ^t'^j ) and by the perturbed geodesic (oc jcj''), following 
from eqn. ilOt and collecting terms quadratic in the gravitational potential, after substitution of eqns. M2>\ and l ll4t : 



3(2) r 



'J 



j(ff{e,x) = -2 r dx' ix-x')- [v.„{x\e,x'\x')^'l]{0,x') + ip,j„,(x\e,x'\x'HlXe^x')\ , (i5) 

Jo X 



where the first-order solutions for the geodesic x^^\9,x') and for the Jacobian matrix j(^^{9,x') need to be substituted. The first term is 
due to lens-lens coupling, as the galaxy image is deformed by the lensing event at x" and is also deformed by further lensing events in the 
large scale structure between and;^'' via ^^1^9, x")- The second term is due to the dropped Born-approximation, i.e. the lensing event at 
x" picks up the derivatives of the potential at the position -I- x^^\9,x') rather than &\.x'9 because of the perturbed geodesic due to the 
integrated lensing effect up to the distance These corrections can be interpreted as nonlocal interactions of the light with the gravitational 
potential and break the symmetry of the Jacobian, because it is no longer simply a projected double derivative of the gravitational potential 
with interchangeable indices. 

Progressing now to the flexion order, the derivative of the Jacobian is expanded to include terms up to third order in<f,S = S"" + 23*'' -l- 
g(2) ^ gO) _|. with the zeroth order solution Sfjl = 0, indicating the absence of lensing effects if the light bundle propagates through a 

constant gravitational potential. The derivative S*'' measures the third derivative of the gravitational potential, 

sfjl{9,x) = -2 ^\x' b(-x')^v,Ax°ie,x'Xx')- (i6) 

The higher order corrections are again derived by substituting the perturbative expansion of all relevant terms in eqn. dlU . and collecting 
terms quadratic and cubic in the gravitational potential. In this way, the second order correction S*^' acquires terms from the perturbed 
geodesic as well as a lens-lens coupling term from the first order Jacobian matrix ,?!*" and a term from the first order Jacobian derivative S*'': 



Bfl{9,x) = -2 r dx'(x-x')- 
Jo X 

The third-order correction S*'' is then given by 



(17) 



<B%{9,x)= - 2£dx'(x-x')— 



'■I' 

Jo 



X 

sX' 



dx' (X -X')^ \<P..n.,^'S^^k + . (18) 

X 

with the necessary subsitutions of the geodesic x'^\9,x'), the Jacobians ^^'\9,x') and the derivatives S^'\9,x'), i = 1, 2. It is not straightfor- 
ward to identify the terms contributing to the correction S*^' as being due to the perturbed geodesic or due to lens-lens coupling, as different 
orders in perturbation theory start mixing. For brevity, we have suppressed the argument (x(9,x'),x') of the gravitational potential if, and the 
argument (9,x') of the geodesic jc*'', the Jacobian matrices ,?l"', ,?l*^' and of S'". 



Born approximation and flexions 5 



deflection angle / 



Jacobian Si = dpi 36 



derivative S = dSildB 



A 



A 



B 



6 O 




Figure 2. Diagramatic representation of the correction terms in the line of sight expressions for the deflection angle fS, the Jacobian J{ and the derivative S. 
The shading of the vertices indicates the order of the derivative of the potential: first order (white), second order (blue), third order (green), fourth order (red). 
The first row assumes that all approximations are in place, the second row describes the Bom-corrections, and the third row sums up the lens-lens couplings, 
all of which are non-local. The number of vertices is equal to the order of perturbation, and therefore to the power of the gravitational potential. 



We will focus on the spectra <S'''S'''> and (S'^'S*^'), and neglect the cross spectrum (S'''©*^^') which is zero for linearly evolving 
scales (being proportional to tp^) and would only be nonvanishing for nonlinearly evolving scales. We focus in this work on the mode- 
coupling generated by lensing itself, and not that introduced by nonlinear structur e formation. The amp litude of the weak lensing bispectrum 
relative to the spectrum is rather small on scales i i 10' (see, e.g. Figs. 4 and 5 in ISchaefer et allbOOSh which leads us to believe that on the 
scales considered here, the (S'''S*~*>-term can be discarded. 

We expect that it is the Bom correction which is more important for the flexion signal than the lens-lens couplings, which can be 
understood from physical arguments: If a light-ray is deflected from its fiducial path by xj,," it encounters changed tidal fields ipijtm introducing 
flexion. The field (pijkm varies rapidly with distance because in general the derivatives of a Gaussian random field fluctuate more rapidly when 
increasing order of the derivative. In contrast, lens-lens terms of the type ^Pimk^^",] and ipi,„S^"\i^ will be smaller because of the line-of-sight 
integration in the computation of and &"\ which averages out fluctuations. For these reasons, we will neglect the cross-correlation 
{S"'S''') as it contains terms of that type. 



3.3 Diagramatic representation of the perturbative expansion 

A visual representation of the perturbative corrections to the lensing effects is given in Fig. |2] This shows graphs corresponding to the cor- 
rections to the lensing deflection angle p (first column), the Jacobian y{ (centre column) and the derivative 2? (right column) for a simplified 
approximative calculation (first row), for the corrections to Bom (second row) and the lens-lens couplings (third row). With all approxi- 
mations in place (first row), the lensing deflection angle p contains gradients of the gravitational potential, the Jacobian ^ contains second 
derivatives and the derivative 2? of the Jacobian contains third derivatives of the gravitational potential, in accordance with Fig.[T] Under these 
approximations, the interaction of the light ray with a gravitational potential is local, and the interchangability of partial derivatives translates 
to the symmetry under index exchange of y{ and S. 

If first-order Bom-corrections to the geodesic are taken into account (second row), the Jacobian y{ is corrected by a term linking the 
third derivative on the fiducial ray with the displacement from the fiducial geodesic (c.f. centre panel in the second row, corresponding to the 
second term in ean.ll5t. and the derivative S measures the contraction of the fourth derivative of the potential along the fiducial ray with the 
displacement (right panel of the second row, depicting the first term of eqn.[T7J. 

Lens-lens coupling (third row), i.e. the interplay between multiple distortions along the line-of-sight, appears in the Jacobian a& & 
contraction between two second derivatives of the gravitational potential (centre panel in the third row of Fig.|2l first term in eqn.|15ll. The 
equivalent corrections for the flexion signal consist of contractions of second and third derivatives of the gravitational potential (right panel 
in the third row): firstly, there is a contraction of the third derivative of the gravitational potential with the integrated second derivative, i.e. 
the Jacobian ^ (second and third term in eqn. I17t and secondly, a contraction of the second derivative with the integrated third derivative 
S (fourth term in eqn. I17t . both of these taking place along the unperturbed ray path. It is important that in this panel both graphs with 
interchanged vertices appear, which is a consequence of the reversability of the light path in geometrical optics. If only one of the graphs was 
present, the reversability would be violated because a specific direction of propagation would be singled out. Due to the fact that both graphs 
with interchanged vertex order are present, no specific direction of propagation is preferred over the other 

So far, we have considered corrections to second order involving products of potentials, which corresponds to the fact that there are two 
vertices in each graph, but with these graphical mles, it is quite easy to generalise the expressions for ^ and S (or even dS/dd) to higher 
orders bearing in mind that we only acquire corrections of order (/>- and therefore would follow the pattern outlined in Fig.|2] 

Quite generally, the lowest order quantities are symmetric under index exchange, reflecting the interchangeability of partial derivatives 
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of the projected gravitational potential, whereas higher-order corrections violate this feature. Hence the Jacobian is symmetric, with 
yi*^' breaking this symmetry. In complete analogy, S*'* is symmetric under index exchange, where again the higher-order corrections S*^* 
and S*^' introduce a symmetry-breaking in the first two indices. Conversely, a useful consistency check is the fact that the expressions 
for S'"* and S'^' derived above respect the symmetry of S under exchange of the last two indices, in accordance with the definition as 
Sijk = dJiijIdOk = dyiikldOj = Bikj- As a last point we mention that the flexion signal has more terms contributing to corrections due to 
dropping the Bom-approximation and lens-lens coupling, compared to the Jacobian. At first order in perturbation theory, there are twice as 
many terms contributing relative to shear, and at second order there are four times as many terms relative to shear. We need to examine in 
later sections whether this increase in number of terms corresponds to an increase in net importance of the corrections at this order. 



3.4 Decomposition in Pauli- and Dirac-matrices 

At this stage, in order to make progress in calculating the corrections, it will be very useful to introduce a decomposition of ^ and S in terms 
of matrices. The Jacobian matrix ^ can be decomposed with the Pauli-matrices cr^, because they constitute a basis for the vector space of 
2x2 matrices, 

3 

If/ = id(2) - ^ - Jl^^'> = ^ a„o-„ = KO-o + y+cTi - ipcTz + j^o-^, (19) 

£1=0 

(c.f. lAbramowitz & Stegunlll972l : lArfken & Webej2005l) . where id(n) denotes the n-dimensional unit matrix. The Pauli-matrices are defined 
as 

/ +1 / +1 \ / -i \ / +1 \ 

= [ +1 )' ^^'=(+1 )' = [ +i J' '^' = [ -1 J' ^20) 

with the properties al = id(2) and tr(tr„) = for o- = 1, 2, 3. Due to the property a-„(rp = id(2)5,^B -I- i^afiyO'y of the Pauli-matrices, the 
coefficients a„ can be recovered by using Oa = ^tr(i/rcr„). 

The coefficients in the decomposition can be identified as the convergence k = oq, the two components of shear 7+ = «[, = and 
the image rotation p = aj. Due to the symmetry of the first order Jacobian matrix gravitational lensing is only able to excite three of 
the four possible image distortions. To first order, it is only possible to observe convergence and the two components of shear. If the Jacobian 
^ beco mes non-symmetric because of the contr ibution then the pa-2 term becomes permissible and the image distortion can include a 
rotation jjain et al.ll2000l : iBacon & Sch"ife3l2009 , for an application). 

In order to carry out an analogous decomposition of the derivative S of the Jacobian, which is a 2 x 2 x 2-tensor, we recast it into a 4 x 4 
symmetric block-diagonal matrix, Sj+2i,k+2i = Sjjt, 

s.-f^r J' (21) 



Si 



where we use the same symbol S as confusion is unlikely. Matrix ( I21b evidently shows the six flexion degrees of freedom, due to the 
symmetry of Sy^.. Sjk can be decomposed with the Dirac-matrices A^: 



(22) 



S = 2 byhy with by = -tr(SAy). 

y=0 

The Dirac-matrices A^ are a generalisation of the Pauli-matrices erg jArfken & Webeil 20051) , of which we select the symmetric subset 



-O-Q / \ -CTl / I -0-3 



(24) 



with the unit element Aq = id(4). These matrices have the properties A^ = id(4), tr(Ay) = for 7 > 1 and form a basis set of the vector 
space of block-diagonal symmetric 4 x 4-matrices, because the erg, /J = 0, 1, 3 are symmetric, cr^ = erg . The fact that there are exactly these 
six Ay-matrices is consistent with the six degrees of freedom of the tensor Sijt, as explained earlier. The selected subset of Dirac matrices 
are a closed group under multiplication, and the product A„A^ of two different A/j-matrices is always traceless such that the decomposition 
by = |tr(SAy) is always possible. 



3.5 Rotations and coordinate-independent representation 

The transformation properties of the S,j-tensor and hence of the coefficients bp under coordinate rotations are complicated, which is due to 
the transformation law 2?,^^ = 2,' 2/ 2t' Rn'Rjj'f^kk'^ii'fk' with three 2x2 rotation matrices This is in contrast to the decomposition of 
the Jacobian into Pauli-matrices, where the coefficient k is invariant under orthogonal transformations of the coordinate system, and the 
two components of shear transform into each other. They can be combined to form the complex shear 7 = 7+ + i7x with the transformation 
property 7 — > yex\>{2iif) under rotations by </). The complexity arises when the tensor S;^^. is mapped onto the block diagonal matrix Sy^ and 
the transformations of the individual indices get mixed. The decomposition of S with Aq and A3 gives rise to a quantity which transforms 
like a vector, i.e. maps onto itself after rotation by 2n, and which is non-zero when the Born approximation is included (see Sect. 4.4). We 
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recognise this as the vectorial spin-1 flexion T\ thus T can be recovered by decomposition of S with Aq and A3. Decomposition of S with 
the remaining A-matrices has a more complex transformation property; however under the Born approximation the decompositions with A2 
and A5 are non-zero, so these are associatied with the flexion Q. The decompositions with Ai and A4 are unactivated in straight light paths, 
and correspond to the twist components. 



4 ANGULAR SPECTRA 
4.1 Angular power spectra 

We now proceed to calculate angular power spectra for the various flexion components, including the perturbation corrections we have 
described above. To lowest order in perturbation theory, the two-point correlation function of the Jacobian matrices and S and hence the 
angular spectra are determined by a projection of the power spectrum of the gravitational potential P^(k). As shown in Sect. [3] dropping the 
Bom-approximation introduces corrections of order 0((p^); these real-space products of (p become convolutions in Fourier-space, so in this 
section we will find terms containing the squared spectrum P^(k) with a mode-coupling kernel. We will truncate the perturbative expansion 
and discard terms of order (p^ because we expect them to be small compared to the corrections 0((f'), and due to reasons of practicality: in 
the power spectra, the corrections 0{tp^) give rise to 6 terms, whereas the corrections 0(ip^) cause a total of 8 + (3) = 64 terms. 
We calculate the correlation function of a„, which are the components of the Jacobian Jl via 

{aje)a;(e')} = ^<tr(,?lWtr„) trmncrpT} = 2 Z (^'■b(f)^h(n} (25) 

a,b iJ 

and apply an analogous relation for the correlation function of the flexion components by, 

{b,({)b',({')) = ^MS({)^y) tr(S(r)A,)*) = ^ Z Z (Soi'(f)^(O) ^.tcAh- (26) 

a.b i.j 

In total, the expression eqn. l l25t contains 4 terms, and eqn. l l26t contains 16 terms, where the signs are given by contraction with cr^i^iCr'i^^ 
and A^i^^A^j., respectively. We introduce the angular power spectra C'^i^-j(C) of the Jacobian 

m„de)^ij(n) = (2nf6o({ - nc-2,j({\ (27) 

and C^i^-j(£) of the derivative S of the Jacobian, 

{S,„(mij(n) = (2Kf6o(e - {')C%(e), (28) 
such that the flexion and twist spectra C^^(C), y, 5 = ... 5, can be defined with (A+ = A, i.e. A' = A and A* = A for the selected subset): 

a.b i.j a,b iJ 

in analogy to the four spectra C'^CO, c,/? = ... 3, describing convergence (a = 0), the two components of shear (a = 1,3) and the image 
rotation (a = 2), 

= i Z Z ^^"-^^^ <^.M^h.- (30) 

a.b iJ 

The power spectrum P^(k, a) of the potential (p = O/c- at the epoch a follows from the comoving Poisson equation A(fi = 3HQCl„J(2a)S and 
is related to the density power spectrum Pg(k, a) by 

__, (31) 

where the Hubble distance = c/Hq makes the fc"'*-factor dimensionless. For simplicity, we only consider linearly growing scales in this 
study, but will make an argument to extrapolate to smaller scales in Sect. 14.41 



4.2 Born-approximation 

In order to project the so urce tenn powe r spectrum Ps(k) in order to estimate the angular power spectrum Ck(£) of the observable, the flat-sky 
Limber-equation is used jLimbeJI 19541) , 

K= r dxW(x)6(x)-^CM= r %W(xfPs(k = e/x), (32) 
Jo Jo X 

with the lensing efficiency function Wix) for sources placed at the comoving distance ^j, 

3Hla„, X, - X D+(x) 

W(x) = —r^ X- (33) 

2c^ Y, a 
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The integration in eqn. l l32t is extended to the comoving distance of the lensed background galaxies Xs- The source term spectrum of the 
shear-related image distortions is given by four-fold differentiation of the potential power spectrum, which results in the angular power 
spectrum C;^,^(^) after Limber-projection, 

= 2' r % f^^l 44 e.tj P,{k = eix,x)- (34) 

Jo ^ \ Xs I 

Spectra of the individual lensing observables can then be obtained by contraction with CT/jjijO'^ji (cf. ean.l27t. This technique is very useful 
for the far larger number of correction terms in the flexion specta. The relevant flexion spectra are derived in complete analogy: the source 
term spectrum follows from six-fold differentiation of the power spectrum of the gravitational potential, with the corresponding implicit 
projection onto Cf^ji^(C), 

Cl2.,.2.,.2,M2.(n = 2' £' ^ (^^) 444 ^,44 P,(k = eix,x)- (35) 

Again, the spectra of the flexion components follow from contraction with Ay_afcAa_,j, according to eqn. l l29t , after recasting the rank-3 tensor 
Sijk into a 4 X 4-matrix Sjt using eqn. ( I2U . 

Due to the transformatio n properties of S, we do not attempt to define a coordinate-i ndependent representation equivalent to the E- 
and B-modes of weak shear ( Stebbins|[T99^ : iKamionkowski et alJI 19981 : IHu & WhitdboOlh . Instead, we provide spectra of the expansion 
coefficients by by choosing the coordinate frame such that the separation vector t coincides with the ^j -axis of the Fourier-frame. In this 
frame, only the component C*(,qq(0 is nonzero, which leads to nonzero projections with Aq, A2, A3 and A5. The vanishing projections are 
those with matrices constructed from cti, which is sensitive to off-diagonal elements in Sat, all of which are zero. This corresponds to the 
absence of a twist mode under the Bom-approximation. In contrast, matrices constructed from ctq and cr^ measure differences between the 
diagonal elements inside each block, as well as between the blocks, and are therefore non- vanishing. 



4.3 Corrections due to the perturbed geodesic, lens-lens coupling and lens-flexion coupling 

In this section, we compute the spectrum of the second order corrections S*^' which are quadratic in the gravitational potential. As 

explained in Sect. |3.2| we expect the terms in to be the dominating contributions to the flexion spectrum. The correction S|J^ to the 

lacobian derivative, at which we truncate the perturbative expansion, will comprise one Bom-term, two lens-terms and one flexion-lens term 
and reads 

^,7*(^'/t)=4 r d/Cr-y)— r '^X" ix' -X'''\vijkm<Pm+y^[vimkV,,ii + Viin,Vml)+\^^ VimVmjX (36) 



after substitution of the first-order expressions for .tJ", andS|]] into eqn. Ml\ . Again, the symmetry 25,^^^ = SjUj remains conserved at all 
orders of perturbation theory. In particular the two lens-lens terms transform into each other under exchange 7 <-> and are both necessary to 
conserve this particular symmetry. 

Consequently, the power spectrum of 2? acquires 6 correctional terms, 3 auto-spectra and 3 cross-spectra. The dependence of the terms 
in S on the square of the gravitational potential maps onto convolutions of the potential in Fourier space, together with contraction of the free 
indices. The d,;^'' -integration will give rise to mode-couplings, weighted by different powers of x' Ix if sach of the source terms S ijk(t,x)'- 

5^(4^) = ^ £^dx'(x-x')vuUnvnAt-t'), (37) 

5^(4^ = J ^ £^ dx' (x-x) j [v,n,k{t')v.,,jit - {') + v,jm{t')v.,,k{{ - n) , (38) 

(4X) = j £ dx'(X-X')i^J^ ViJOVrndt-n 



s 



(39) 



These source terms can be combined to yield the spectrum: 

(S„tA{,x)S:Af',x')) = ^^Soix - x')Sd({ - nP.bciik({,xl (40) 
X 

which can then be projected for obtaining the angular spectra: 

Ct^2..c^2..jU,M2,(t) = 4- ^ (^y^) P'"'c,jk{l,Xy (41) 

The resulting power spectra between the Born-, lens- and flexion-terms are obtained by applying the Wick theorem (cf . lDurreill2008l) to the 
4-point correlation of the gravitational potential. It is worth noting that the various resulting spectra have very similar forms, differing only 
in the extent of mode coupling and dependence on £: 

CU^'X) = J^[e'(e-nf{VXf:f'A)Mo(ej',x\ (42) 



Born approximation and flexions 9 




Figure 3. Mode coupling functions M^(C, {',x), for 1 = ... 4 (from top to bottom) used in the corrections to the angular flexion spectra, computed for a 
distance of ^ = 1 Gpc//i, which corresponds to a redshift ofz- 0.36 in ACDM with £2„, = 0.25. The angle <f>(> between t and t' is fixed at a value of 7t/2. 

Pl%k(f'X) = J ^ [f'(f - nf [Qt - t'Ut - t'\ l\{t - t')j{l - f ),) M,{t, t\x\ (44) 
The corresponding cross spectra, where a symmetrisation has been applied, are given by 

WAH - OAk + e'aVAAjd - Ok) M^e, t',x\ (45) 

p'TcjS^^x) = J ^ [t'(e - nf [Q{ - nAAd - njit - Ok + w - OAi^ - OjV - (')k+ 

- ndf - £')A(£ - e')A + i'aV - e'ue - OAAM - ^')*) ^3^^' m 
pHUAx) = / ^ [^'(^ - [VAc m - Ojie - e')t + r„(e - rw - He tAA^ e',x)- (47) 

Common to all spectra is the mode coupling function Mg, which arises from the fact that the corrections are proportional to the square of 
the gravitational potential in real space, which translates into a convolution in Fourier-space. The d-T -integration is carried out in polar 
coordinates, (ft' = E'AE'A^p with t aligned along the f^-axis such that t = ^(1, 0) and f = ('(cos (pf, sirn^f). The mode coupling Mg({, t ,x) 
reads: 

M,(e,e',x) = 2£ ^(j(-x'f{jj PAf'\A)PAf-f'\'X')- (48) 

The mode coupling kernels Mg(£,{',x) are depicted for a fixed distance ;f = 1 Gpc//j and a fixed angle ipc between i and {' of 7i/2 in 
Fig-El for the range of indices q = . . .4 considered here. The strongest mode coupling is between adjacent f-modes and drops rapidly with 
increasing distance in Fourier-space. Fig.[3]confirms our qualitative argumentation at the end of Sect. |3.2| concerning the relative magnitude 
of the Bom- and lens-lens terms: M^(f, f',x) assumes the largest values for small q, and in particular for ^ = 0, which is the spectrum of the 
Born-correction only, and decreases if lens-lens terms are considered instead. This provides a quantitative argument for our expectation that 
the largest contributions originate from the Bom-terms in the (S*^'S'^'>-spectrum. 



4.4 Corrected angular flexion spectra 

The spectra of the flexion coefliicients are shown in Figs.|4]and[5] where the lensed background galaxies are assumed to reside at a redshift 
of Zs = 0.9 (~ 2.2 Gpc//? in ACDM), which corresponds to the median redshift of the EUCLID galaxy sample. The corrections to the 
spectra comprise contributions due to the dropped Born-approximation, lens-lens and lens-flexion couplings, to linear order in the perturbed 
geodesic. If all approximations are in place, the spectra involving decompositions with Ai and A4 are equal to zero, and all other spectra 
are identical. Corrections arising if the weak flexion signal is evaluated along the true geodesic instead of the Bom-approximated geodesic 
amount to 10"'' on the smallest angular scales considered, with the difference increasing towards larger angular scales. Nonlinear corrections 
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Figure 4. Flexion angular power spectra Cys{t) assuming the Bom-approximation (thick solid line) and corrections with a perturbed geodesic (thin solid lines) 
for y, (5 e |0, 2, 3, 5|. The background galaxies are all placed at a redshift of Zs = 0.9, and the cosmic variance eiTor ACyi(^) is plotted for comparison (thick 
dashed line). 

to the CDM spectrum P{k) on small spatial scales would increase the amplitude of the power by approximately 1.5 dex jSmith et al.ll2003h . 
and increase the lensing correction amplitude by 3 dex at most (since the corrections are proportional to the square of the power). Hence the 
relative correction is increased from 10 ■* (in the linear regime) to 3 x 10"^ in the very nonlinear regime; thus we conclude that for cosmic 
flexions, the Born- approximation is fulfilled to a very high degree of accuracy. This is despite the fact that there are more contributing terms 
in the flexion correction in comparison to shear co rrection; i.e. we find that flexion correction te rms are proportionately somewhat smaller 
than the shear corrections terms already calculated IShapiro & Cooravll2006l ; ICoorav & Hij|2002l) . 

Note that on very small scales one starts to probe the internal structure of haloes instead of the smooth Gaussian fluctuations of the 
cosmic density field. In this case, it would be more appropriate to replace the C DM spectrum P(k) with a description following the halo- 
model (for a review on the halo model including its application to weak lensing, see ICoorav & Sheflj2002l) . The more accessible and possibly 
more relevant quantity would be the 1-halo term, for which one would need to compute the above discussed correction in the case of lensing 
on a single, spherically symmetric matter distribution. The lens-lens couplings, however, would be difficult to incorporate into this model due 
to the correlatedness of the individual deflecting haloes along the line of sight. 

Comparing the magnitude of the corrections to the uncertainty due to cosmic variance, 

AC,.(0=^|I|ITc,,(,), (49) 

one sees that, for a Euclid-like survey with f^^y = 1/2, the corrections are much smaller than cosmic variance for the range of multipoles 
considered here. Note that the newly excited twist degrees of freedom associated with the decomposition with Ai and A4, reach amplitudes in 
their power spectra which are slightly smaller than the corrections to flexion power spectra. Hence twist from a cosmic origin will be difficult 
to detect; if it is present in surveys, it is likely that an origin from systematic effects will dominate. 



5 CONCLUSION 

This paper treats corrections to the angular spectra of weak cosmic flexions due to nonlocal interactions of the light ray with the cosmic large 
scale structure, which arise if the Bom-approximation is dropped and if lens-lens and lens-flexion couplings are active. We have derived 
corrections to the spectra by perturbing the geodesic to first order. All integrals needed for computing the geodesic corrections can be 
represented in a graphical way, which we have summarised in Figure [2l For illustration, we have carried out the calculations for a ACDM 
cosmology for the planned EUCLID weak lensing survey, which reaches a median redshift of 0.9. 

The corrections violate the symmetry of the Jacobian and the derivative S due to their nonlocal nature. Concerning the flexion signal, 
these fall into three categories: the perturbed geodesic makes the light ray experience gradients in the gravitational potentials at a different 
position relative to a fiducial straight ray, and solving the implicit propagation equations for ^ and S with a perturbation series give rise to 
lens-lens couplings and lens-flexion couplings, where the integrated lensing signal is coupled to a derivative of the gravitational potential. 

We have explained how the perturbative expansion of the line-of-sight integrals is a perturbation series in the gravitational potential. 
Powers in the gravitational potential give rise to mode-couplings in Fourier-space when computing angular spectra. Naturally, the mode 
coupling is strongest between adjacent Fourier-modes and drops rapidly with increasing distance in Fourier space. 

We have decomposed the spectra in terms of Pauli- and Dirac-matrices, respectively, which allows us to separate the newly arising 
distortion modes, and to identify the vectorial flexion. We have derived the perturbative corrections to all spectra of the flexion components, 
and have shown that they are smaller than the Bom-approximated spectmm by approximately four orders of magnitude, well below the 
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Figure 5. Angular flexion power spectra Cys(f) involving y, i? e {1,4) (thin solid lines), which vanish if the Born-appoximation is applied, but are non-zero 
from our calculation using a perturbed geodesic. We plot for comparison the flexion spectrum Coo({) (thick solid line) with its cosmic variance eiTor ACoo(^) 
(thick dashed line) assuming the Born-approximation (compare Fig.|4). The background galaxies are assumed to reside at a redshift of Zs = 0.9. 



cosmic variance limit. The spectra of the newly arising twist components are of similar magnitude compared to the corrections, albeit 
slightly weaker. All corrections are of similar relative size compared to analogous corrections to weak cosmic shear spectra. 

For the case of weak cosmic flexions, we conclude that the Bom-approximation remains an excellent approximation. We expect that 
geodesic corrections are larger in the case of cluster lensing, as the second derivatives of a cluster lensing potential are stronger than those 
arising in the large-scale structure. It remains to quantify the importance of the flexion corrections in that context, together with the amplitude 
of the newly excited twist-degrees of freedom associated with Ai and A4. 
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